Spin and Orbital Order of the Vanadium Spinel MgV 2 04 
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We present a unique study of the frustrated spinel MgV^CU which possesses highly coupled spin, 
lattice and orbital degrees of freedom. Using large single-crystal and powder samples, we find a 
distortion from spinel at room temperature (space group F43m) which allows for a greater trigonal 
distortion of the VOs octahedra and a low temperature space group (I4m2) that maintains the mirror 
plane symmetry. The magnetic structure that develops below 42 K consists of antiferromagnetic 
chains with a strongly reduced moment while inelastic neutron scattering reveals one-dimensional 
behavior and a single band of excitations. The implications of these results are discussed in terms of 
various orbital ordering scenarios. We conclude that although spin-orbit coupling must be significant 
to maintain the mirror plane symmetry, the trigonal distortion is large enough to mix the 3d levels 
leading to a wave function of mixed real and complex orbitals. 



PACS numbers: 75. 25. +z, 75. 10. Jm, 61.12.Ex 

Geometrically frustrated magnets are characterized by 
competing interactions resulting in a highly degenerate 
lowest energy manifold. In many cases the degeneracy is 
eventually lifted at low temperatures by a lattice dis- 
tortion. In compounds where the magnetic ions also 
possess orbital degeneracy, orbital-ordering can influence 
the exchange interactions and lift the frustration. The 
vanadium spinels AV2O4, where A is diamagnetic Cd 2+ , 
Zn 2+ , or Mg 2+ provide ideal systems to study the inter- 
actions between spin, lattice and orbital degrees of free- 
dom pp. In these compounds the magnetic V 3+ -ions pos- 
sess orbital degeneracy and form a frustrated pyrochlore 
lattice with direct exchange interactions between nearest 
neighbors providing a direct coupling to the orbital con- 
figuration. The interplay of orbital and spin physics has 
been studied in other systems like the perovskite vana- 
dates (RVO3, R is a rare earth) which show a strong 
correlation between orbital ordering and magnetic struc- 
ture [2]. However in these compounds the couplings are 
unfrustrated and indirect, occurring via super exchange 
through oxygen. In contrast the magnetic structure and 
excitations in the spinel vanadates are more sensitive to 
orbital ordering and thus characteristic of it. Further- 
more the additional component of frustration allows for 
the possibility of exotic ground states. Indeed the nature 
of the ground state in the spinel compounds has gener- 
ated intense theoretical interest over the past eight years 
(3HZ] but has remained an unresolved experimental issue 
which we will address in this paper. 

The electronic configuration of V 3+ is 3d 2 leading to 
a single-ion spin 5=1. Each V 3+ -ion is located at the 
center of edge sharing VOq octahedra which create a 
crystal-field that splits the d-orbitals and lowers the en- 
ergy of the three t2 g orbitals by approximately 2.5 eV [3]. 



These levels are usually assumed to be degenerate (ignor- 
ing the small trigonal distortion which will be discussed 
later) so that the two (i-electrons of V 3+ randomly oc- 
cupy the three ti g orbitals, introducing an orbital degree 
of freedom. This class of compounds undergoes a struc- 
tural transition (Tg) from cubic to tetragonal which par- 
tially lifts the orbital degeneracy. A magnetic transition 
to long-range antiferromagnetic order occurs at a lower 
temperature (Tn). This paper describes the first com- 
prehensive single-crystal investigation with supporting 
powder measurements of the crystal structure, magnetic 
structure and magnetic excitations of MgV204 which al- 
lows the orbital ordering to be determined and reveals 
the interplay between frustrated magnetism and orbital 
physics. 

All measurements of the AV2O4 series have so far 
been performed on powder samples. Z11V2O4 is the 
most extensively studied member of the series. Pow- 
der neutron diffraction measurements reveal a cubic to 
tetragonal structural phase transition at Tg=51K where 
the tetragonal c-axis becomes suppressed and the space 
group is 74 1 / amd [8] . Long range antiferromagnetic or- 
der occurs at Tn=40K with an ordering wave vector of 
k=(0,0,l). Powder inelastic neutron scattering (INS) re- 
veals a change in line shape of the magnetic fluctuations 
at Tg attributed to a change from three-dimensional to 
one-dimensional coupling [9 . The tetragonal distortion 
at Tg causes the t 2g levels to split into a lower d xy level 
and a twofold degenerate d yz /d xz level as shown in Fig. [I] 
One of the two d-electrons occupies the d xy level lead- 
ing to antiferromagnetic chains in the [l,l,0] cu bic an d 
[1,1,0] cu bic directions. The second electron occupies ei- 
ther d yz or d xz or some combination and is responsible for 
the inter-chain coupling along [l,0,l] cu bi c and [0,l,l] cu bic- 
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FIG. 1: VC>6 octahedra along the chain direction in both the 
cubic and tetragonal phases (distortions have been exagger- 
ated). The trigonal distortion slightly splits the £29 levels in 
the cubic phase. In the tetragonal phase the oxygen octahedra 
are compressed causing large splitting of the t2 a levels. 



Models proposed for Z11V2O4 either place the second elec- 
tron in a real combination of the d-orbitals thus maintain- 
ing the orbital quenching or in a complex combination re- 
sulting in unquenched orbital moment. We will consider 
the three main models. 

The first model is a real orbital ordered model (ROO) 
where the Kugel-Khomskii exchange Hamiltonian is in- 
troduced on the pyrochlore lattice. The ground state con- 
sists of alternating occupation of the d yz or d xz orbital in 
successive layers along the c-axis [4] . It gives rise to weak 
and frustrated ferromagnetic inter-chain coupling and no 
spin anisotropy. However, it implies that the space group 
is J4i/a, which has lower symmetry than that of ZnV204. 
Therefore a complex orbital order model (COO) was pro- 
posed that is compatible with the space group IAi/amd 
of ZnV204 which has additional mirror and diamond 
glide planes. In this model the second electron occu- 
pies the orbital (d yz ± id xz ) which has unquenched or- 
bital moment of L=l [5]. The resulting inter-chain cou- 
pling is frustrated and antiferromagnetic and is expected 
to be larger than in the ROO case. It is favored when 
the spin-orbit coupling is strong in comparison to inter- 
chain exchange interaction because the crystal can lower 
its energy by aligning L opposite to S. This scenario 
satisfies the spin-orbit interaction resulting in a reduced 
moment of 1 /is compared to the spin only value of 2 
in ROO. In contrast the ROO model is favored when the 
inter-chain exchange interaction is dominant because it 
reduces the electron overlap and thus the energy of the 
frustrated bond. The third model considers the prox- 
imity of ZnV204 to the itinerant-electron boundary [7]. 
Electron band structure calculations imply a structural 
distortion where the vanadium- vanadium bond distances, 
d(V-V), alternate along the inter-chain directions. The 



shorter bond is predicted to be 2.92 A , below the limit 
for itinerancy. It is therefore described as a homopolar 
bond (HPB) characterized by partial electronic dereal- 
ization which results in ferromagnetic alignment of spins 
and release of inter-chain frustration. In this model the 
wave function would be the real combination of orbitals 
{dy Z ± d xz ) with no net orbital moment. 

Further distortions of the V06 octahedra besides 
the tetragonal distortion considered by these models 
are possible in the AV2O4 compounds. These could 
mix the d-orbitals and completely remove their de- 
generacies resulting in a general wave function e.g. 
ad X y+/3d xz +jdyz+e(eg-orbitals). Indeed a trigonal dis- 
tortion is already present in the cubic phase of ZnV2 04. 
If the resulting splitting dominates over the exchange in- 
teraction or spin-orbit coupling then a wave function sim- 
ilar to that calculated for the related material MnV204 
may result where a, /3, 7 and e are all non-zero and 
mostly real [10] . 

The lack of large single crystals has prevented the full 
determination of the structural and magnetic behavior of 
the AV2O4 series. Here we report the first single- crystal 
study of MgV2 04. This material was first investigated by 
Plunder and Tardieu who observed magnetic order at low 
temperatures [11] . Magnetic susceptibility, specific heat 
and powder x-ray diffraction measurements by Mamiya 
et al. showed two transitions at Tg=65K and Tn=42K 
with the system changing from cubic to tetragonal below 
Ts |12j . Here we describe x-ray and neutron diffraction 
and INS in the cubic, tetragonal and magnetic phases. 
The data is used to distinguish between the theoretical 
pictures and is found to favor a mixture of real and com- 
plex orbitals. 

To precisely determine the crystal and magnetic struc- 
ture of MgV2 04 we carried out two complementary 
diffraction experiments. The 4-circle diffractometer E5 
at the BER II reactor (Helmholtz Zentrum Berlin) was 
used with a wavelength A=0.89A along with a 3 mm 
Er filter and A=2.38A with a PG filter to measure a 
large (187 mm 3 ) high-quality single-domain crystal [T5] . 
Powder x-ray measurements were made using the ID31 
diffractometer at the European Synchrotron Radiation 
Facility with A=0.5A at 10 K and A=0.4A at 300 K. 
Refinements were performed using the FULLPROF [T3] 
and Xtal programs [TS] . 

First we consider the room temperature results where 
the system has cubic symmetry. The single-crystal neu- 
tron diffraction revealed the presence of low intensity re- 
flections of the type (hOO) where h=2n and (hkO) where 
h, k=2n. These reflections are forbidden in the spinel 
and reveal the loss of the diamond glide plane (Fig.|2|. 
They indicate that the space group is F43m. The trigo- 
nal symmetry of the vanadium site characteristic of the 
spinel is maintained. Refinement of the room tempera- 
ture x-ray and neutron diffraction yield the results listed 
in tableU The vanadium ions were found to remain at the 
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TABLE I: Lattice parameters and fractional coordinates in 
the tetragonal (6 K) and cubic (300 K) phases. In the cubic 
phase Mg sits at 4a(0,0,0) and 4c(|,|,|). Both O and the 
V sites are 16e(x,x,x). In the tetragonal phase Mg sits at 
2a(0,0,0) and 2c(0,|,») while the O and V sites are 8i(x, 0, z). 
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FIG. 2: Neutron single-crystal diffraction results. Intensity 
of both the main Bragg reflections, (222), and the super cell 
reflections, (002), show a jump at Ts- At Tn, magnetic Bragg 
reflections at structural absences such as the (001) appear. 



ideal spinel position. Only the oxygens distort the lattice 
from spinel to increase the trigonal distortion around the 
vanadium ion. 

Low temperature .T-ray powder diffraction showed that 
at Tg no new Bragg reflections appear but that peaks 
split with an intensity ratio of 2:1 indicating a cubic to 
tetragonal structural phase transition where the c-latticc 
parameter is suppressed c/a=0. 9941(1). This value is 
in good agreement with c/a=0.993 found previously [12]. 
For the single-crystal neutron diffraction we ensured a 
single tetragonal domain by applying a uniaxial pressure 
of approximately 0.5 MPa along [001]. Fig. [2] shows a 
jump in the intensity of Bragg reflections at Ts=64K 
due to a change in extinction parameters. The low tem- 
perature data set matched the space group IAm2. This 
space group is characterized by a mirror plane and two 
fold rotational axis which ensures that the spatial wave 
functions of the electrons consists of equal amplitudes of 
d xz and d yz , a combination favored by spin-orbit cou- 
pling. Again we found the vanadium was not displaced 
from ideal spinel position. Therefore d(V-V) is not sig- 
nificantly affected by the tetragonal distortion. Within 
the tetragonal plane there are two inequivalent intra- 
chain bonds both of distance d(V-V)=2. 9802 Awhile be- 
tween the planes there are four inter-chain bonds of 
d(V-V)=2.9714A. This result discounts the HPB model 
for MgV2 04 which is characterized by alternating bond 
lengths, the smaller of which should be less than the limit 
for itinerancy of 2.94 A. 

New Bragg reflections appear in the neutron diffraction 



below Tn=42K and can be indexed with the ordering 
wave vector k=( 0,0,1) indicating an antiferromagnetic 
spin alignment along the chains as found previously |11) . 
The spins point predominantly along the c-axis with a 
moment of fi z =0A7(l) fi B - A tilt of 8(1)° from the c- 
axis, similar to the tilt of the oxygen octahedra along the 
chains, was identified by a finite intensity on the (001) 
reflection (Fig[3]a). The moment is greatly reduced from 
the 2 /xb expected for ROO in the case of full magnetic 
ordering and even reduced from the 1 /xb expected for 
COO. It is smaller than in CdV204 or ZnV204 which 
have moments of 1.19/ib and 0.63 /iB, respectively. 

We made INS measurements on the IN20 spectrom- 
eter with the Flatcone detector at the Institute Laue 
Langevin using the detwinned crystal at 2K with the 
tetragonal plane in the horizontal scattering plane. We 
used a Si monochromator and analyzer with a fixed final 
wave vector k^=3A -1 (Ej=18.6 meV). At low energies 
we observe intense scattering along two sets of stripes 
in the tetragonal plane which are perpendicular to each 
other (Fig.[3]b). At higher energy these stripes split into 
two branches. The stripes indicate that modes disperse 
strongly along the straight weakly-coupled chains of an- 
tiferromagnetic spins confirming occupation of the d xy 
orbital. Further measurements at the MERLIN time-of- 
flight spectrometer at ISIS using a powder sample and 
incident energy of 200 meV reveal a single band of exci- 
tations extending up to 65 meV. No magnetic scattering 
was observed between 65meV and 110 meV as shown in 
Fig.|3]c. The magnetic excitation spectrum has been cal- 
culated by Perkins et al. for both the COO and ROO [BJ. 
In the COO there are acoustic and optic bands. The optic 
band is predicted to have possess comparable scattering 
weight to the acoustic band and is a direct consequence 
of unquenched orbital moment. It is expected to lie at a 
higher energy, well separated from the acoustic band due 
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FIG. 3: a) The tetragonal unit cell with only vanadium (blue) and oxygen (red) ions shown. Light/dark blue arrows indicate 
the magnetic moment directions of the two perpendicular antiferromagnetic chains, b) Single crystal INS data measured at 
lOmeV and 1.5 K labeled in tetragonal notation. The intense stripes indicate highly one-dimensional perpendicular chains of 
spins, c) Plot of INS with background subtraction from a powder sample summed over two |Q| ranges. Magnetic signal appears 
in the low |Q| range while at higher |Q| magnetic intensity is suppressed due to the magnetic form factor of the V 3+ -ion. No 
magnetic scattering is observed above 65meV. Lines indicate the bandwidths predicted for the COO model and ROO model [BJ. 



to the spin-orbit coupling. For the ROO model gapless 
acoustic spin-wave modes are predicted which disperse 
up to 40 meV along the chain but are only weakly dis- 
persive between chains. The magnetic excitation spec- 
trum of MgV204 appears to be in better agreement with 
the ROO model. We do indeed observe highly dispersive 
excitations which is contrary to expectations for largely 
local orbital excitations in the COO model. 

These results have implications for the orbital and 
magnetic behavior of MgV2 04 which we now discuss. 
The HPB model is discounted since we have no alter- 
nation of the inter-chain d(V-V). The COO model has 
unquenched orbital angular momentum that reduces the 
total magnetic moment [5] . We do indeed observe a 
strongly reduced ordered moment which is even smaller 
than the 1 hb expected. The 8° tilt of the moment im- 
plies some anisotropy indicating the presence of spin- 
orbit coupling. However, the measured magnetic exci- 
tations are not consistent with a large unquenched or- 
bital moment but are in better agreement with the ROO 
model. Here the orbital moment is quenched and the to- 
tal moment should be 2 /ib however it can be strongly 
reduced due to quantum fluctuations known to suppress 
ordering in one-dimensional systems with small and frus- 
trated inter-chain coupling [HE)]. ROO was the original 
model proposed for AV2O4 systems. It was discounted as 
a description of ZnV2 04 because the mirror plane per- 
pendicular to the tetragonal a or 6-axis constrains the 
occupation of the d yz and d xz orbitals to be equivalent 
on every site. This constraint also holds for MgV204. 

In conclusion, the electronic wave function of the 
V 3+ -ions in MgV204 cannot simply be composed of 
purely real or complex orbitals. Consider a general 
wave function for the AV2O4 compounds of the form 
ad xy +Bd xz +^e l ^d yz . When the tetragonal distortion is 
much larger than the trigonal distortion, which tends to 
split and mix the t2 S -orbitals, we expect j3, 7<Cl 



for the first electron and a<§;l, /3,7~l/v2 for the sec- 
ond electron. In the case of ZnV2 04, density function 
theory (DFT) calculations within the IAi/amd space 
group show the trigonal distortion is small and spin- 
orbit coupling dominates giving rise to COO (i.e. 0=tt/2, 
/3=7=l/\/2 for the second electron) [T7]. In contrast for 
MnV204, where the space group is I4±/a and both the 
diamond glide plane and mirror plane are lost, the trigo- 
nal distortion is large. As a result both electrons occupy 
a mostly real mixture of d-orbitals, </>«0 and B^i |10j. 
The influence of spin-orbit coupling is weak in compari- 
son and the orbital moment is small at 0.34 /ie [l0|, I18j. 
We suggest that the wave function of MgV2 04 sits be- 
tween these two examples. The octahedra around the 
V 3+ in MgV2 04 are significantly distorted which would 
probably split and mix d-levels. The influence of the spin- 
orbit coupling is substantial because the mirror plane is 
maintained constraining 8—j nevertheless the large ex- 
citation band width indicates that the orbital moment is 
not maximal, i.e. /3<l/v2 and 4>^tt/2. We suggest that 
the trigonal distortion also plays a significant role in de- 
termining the wave function and that DFT calculations 
which consider the full VOg distortion would indicate a 
mixture of real and complex orbitals. 
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